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NOTE ON TEE AERODYNAMIC HEATING OF AN OSCILLATING SURFACE 

04 By Simon Ostrach 
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SUMMARY 

An analysis of the temperature distributions in a fluid over an 
osc ill ating surface with heat transfer is made and associated heat- 
transfer parameters are compared with those for the case of conduction 
at a stationary surface with the same initial temperature potential. 

It is found that the heat tr ans fer for the oscillating surface can he 
considerably different from that for conduction alone. The effect of 
the surface oscillations on the thermal state of the fluid is studied 
by means of average static- or total -temperature defects, and it is 
demonstrated that the oscillations could alter the fluid temperature 
* appreciably. 


INTRODUCTION 

The increased improvement of present-day propulsion systems and 
the development of new propuls ion systems have posed numerous new prob- 
lems in the field of heat transfer. Elucidation of unusually high heat- 
transfer coefficients which are apparently encountered in unsteady flows 
and means of increasing heat-transfer coefficients under given condi- 
tions are greatly desired. As a preliminary attempt to gain insight 
into such problems, it seems worth while to consider the heat-transfer 
aspects of the classical problem wherein the fluid motion Is induced by 
oscillating a conducting surface axially in viscous fluid. In particu- 
lar, the effect of disturbing the equilibrium (steady state) conditions 
after the periodic motion of the fluid has been established will be 
studied. The temperature distributions in the fluid are determined as 
exact closed-form solutions of the energy equation pertinent to the 
problem and, hence, related heat-transfer parameters can be compared 
with those for a stationary surface to demonstrate the effect of the 
surface oscillations. Other exact solutions of the energy equation for 
somewhat analogous problems are presented in references 1 to 3. 
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ANALYSIS 


The laminar motion of an incompressible viscous fluid induced by 
the axial (longitudinal) oscillations of a bounding plane surface doubly 
infinite in extent is veil known and is reported in numerous sources 
(see Schlichting, ref. 4, for example). The equations of motion for 
this case, assuming constant physical properties, reduce to the class- 
ical heat -conduct ion equation 


du B 2 u 


CD 


where u is the velocity component parallel to the surface, t denotes 
the time, v is the kinematic viscosity, and y is the coordinate nor- 
mal to the surface. (All symbols are defined in the appendix.) The 
associated boundary conditions are 

u(0,t) » U cos nt (2) 

and 

u(-,t) » 0 (3) 

where U denotes the amplitude and n the frequency of the oscilla- 
tions . The velocity distribution, after sufficient time has elapsed for 
the periodic motion to be established, is given by 


u(y,t) = U exp -y/n/2v yj cos ^nt - -y/n/Zv 3^ 


The appropriate energy equation is 


dr 

^t 


d^T v f ckA 


( 4 ) 


( 5 ) 


where T denotes static temperature, a is the thermal diffusivity, and 
Cp is the specific heat at constant pressure. Note that the last term 

on the right of equation (5) is due to aerodynamic heating and describes 
the increase in fluid temperature caused by the oscillations. In order 
for the problem to be physically reasonable, however, after a long per- 
iod of time the system should be in a state of thermal equilibrium (that 
is, in essentially a steady state) and it is thus necessary to cool the 
surface. The thermal equilibrium state can be described by first solv- 
ing equation (5) subject to the boundary conditions 
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and 


T(-,t) = T m 


(?) 


that is, the osc ill at ing surface Is maintained at some uniform tempera- 
ture T^.. It is unnecessary here to specify any initial conditions, for 
they w ill not influence the results after a sufficiently long time. 


After the explicit form of the nonhomogeneous or aerodynamic heat- 
ing term in equation (5) is determined hy means of equation (4), a par- 
ticular solution of equation (5) was found to he 




U^r 

4C P 


exp 


( — / pyi Ar *tr ] 

, cos(2nt - -\/2n/v y) 

^ ^^n/v yj 

1 + 2 - Pr 


( 8 ) 


where Pr =* v/ot, is the Prandtl number. The complementary function to 
be added to equation (8} in order that the complete solution satisfies 
the boundary conditions given by equations (6) and (7) is 


T c s* T, 


W + ^ £ + ^ c ^(2nt - y) 


erf(-^| 

\^ 2 -Jafcj 


0) 


where erf( ) is the error function and is defined as 

. £ 

erf (p) = 




exp(-r^) dr 


so that 

erf(o) = o 

and 

erf(") =* 1 

The complementary function is synthesized from well-known solutions of 
the heat equationj for example, the term containing the product of the 
exponential and trigonometric functions corresponds to the solution of 
the heat -conduct ion equation for the case where the temperature on a 
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stationary surface is specified to "be a harmonic function of time (see 
ref. 5, p. 47), and the two error-function terms combined correspond to 
the solution of the problem of unsteady heat conduction to or from a 
stationary surface due to a temperature difference or potential of 
U 2 Pr - 

T„ - 4 c“ “ ® w - Note that the error -function terms describe the tran- 
sient heat transfer which will vanish for large time. The temperature 
distribution in the equilibrium state (that is, after a long time) should 
be independent of time except for periodic terms and is given by the sum 
of equations (6) and ( 9 ) with the error-function terms vanishing as 



The relation for thermal equilibrium between the surface temperature T w 
and the ambient temperature T„ can be found by evaluating equation 
(10) as y • and is given by 



( 11 ) 


Thus, as was anticipated, it can be seen from equation (ll) that for 
thermal equilibrium the surface must be cooled to compensate for the 
aerodynamic heating, and the temperature distribution above a surface 
oscillating for a long time in a fluid at temperature T^ is given by 



Now that the fully developed state of the fluid is completely de- 
scribed (by eqs. (4) and (12)), consideration can be given to the prob- 
lem at hand: namely, the effect of the oscillations on the heat trans- 

fer if, at seme time (say t = 0) after the periodic motion (steady 
state) is established, the surface temperature is fixed at some temper- 
ature T' w different from that given by equation (ll) . The problem to be 
solved to answer this question is mathematically identical to the one 
for thermal equilibrium except that T w replaces T w in the boundary 
conditions and that the solution of the present problem must match that 


f 


3169 



MCA TN 3146 


5 


given by equation (l2) for t = 0; that is, the proper initial condi- 
tion must be satisfied by the temperature distribution in order to dem- 
onstrate the transient effects properly. It was previously shown t hat 
the sum of equations (8) and (9) 



satisfied the differential equation (eq. (5)) and the bo undar y condi- 
tions (eqs. (6) and (7)), where T w is now replaced by T w . Evalua- 
tion of equations (12) and (l3) at t = 0 shows that they are identi- 
cal so that the proper initial condition is satisfied and thus equation 
(13) represents the temperature distribution over an oscillating sur- 
face which is maintained at a constant temperature T w< Because T w 
is a surface temperature different from that for thermal equilibrium, 
the solution given by equation (l3) is of greatest physical significance 
only as long as the transient terms are Important. This of course, is 
the problem of primary interest herein. 

Although the solution given by equation (13) was developed for an 
incompressible viscous fluid with constant property values, it is equally 
valid in the case of a compressible viscous fluid if the boundary-layer 
assumptions are made, if the Prandtl number and the product of p and 
p are taken to be constants, and if y is replaced by q where 

T| = r <3y 
'-'0 

This follows because the von Mises transformation (see ref. 6, for 
example) under these assumptions reduces the compressible boundary- 
layer equations to the forms of equations (l) and (5). 

Any pertinent heat-transfer quantities can be obtained from equa- 
tion (l3) . However, since the time dependence itself is usually not of 
primary practical importance, the temperature gradient at the wall 
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(which is a measure of the heat transfer) will be averaged over a period. 
This average gradient is given by 



where tg indicates the beginning of any cycle . For the equilibrium 
state, an analogous equation is obtained from equation (l2): 



(15) 


This average temperature gradient is due entirely to the aerodynamic 
heating and is also essentially (that is, except for the multiplicative 
constant k, the thermal conductivity coefficient) equal to the work 
done per cycle in oscillating the surface. The latter result can be 
independently verified by computing the work done per cycle from the 
shear stress on the plate. 


RESULTS AND DISCUSSION 

Now that the temperature distribution and its gradient at the sur- 
face are known, the effect of the surface oscillations on the heat 
transfer can be obtained by comparing any appropriate parameters with 
the corresponding ones for the case of pure conduction to or from a 
stationary surface subject to the same initial temperature difference 
(T„ - T w ) . The specific formulation of the parameters depends on the 

particular configuration considered. With this in mind, the parameters 
considered herein will fall into two categories: total parameters, that 

is, those pertaining to configurations in which the oscillations are 
inherent} and net parameters for configurations in which the work done 
in oscillating the surface is assessed to the system. For example, the 
average total heat-transfer coefficient as given essentially by equation 
(14) can be compared with that for conduction to or from a stationary 
surface at an initial temperature difference (T„ - Ty.) . The temperature 
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distribution for the latter case is given by the sum of the first and 
last terms of equation (l3) (see ref. 5, p. 34) and the associated aver- 
age temperature gradient is given by merely the last term, of equation 
(l4) . Therefore, the ratio of these average total heat -transfer coef- 
ficients R^ is given by 



It can thus be seen that the osci ll ations could lead to higher average 
total heat -transfer coefficients depending on the sign of (T„ - T ) and 
on the relative magnitude of the two terms in the bracketed part of equa- 
tion (16); for example, in the limiting case just after the temperature 
T w is imposed on the surface, that is, tQ -*■ 0, equation (l6) reduces to 

E io - 1 - - *> W 

so that for the case where the ambient temperature is greater than the 
surface temperature (T„ >T w ), and Pr < jt , then R 10 >1; that is, the 
oscillations increase the average total heat-transfer coefficient. Con- 
versely if T w > T„., R-jj-j would be larger than unity if Pr > it 2 . 

It should be reiterated that the parameters for comparison depend 
upon the configuration considered. Therefore, if in the specific appli- 
cation the net energy obtained from this system (rather than the total 
as before) is of interest then the total heat transfer must be properly- 
assessed to take account of the work done in oscillating the surface. 
Hence, the net heat transfer for this case (as represented by the dif- 
ference between eqs. (l4) and (lb)) will be compared with that for con- 
duction to or from a stationary surface. The ratio Rg of the average 
net heat transfer for the oscillating surface to that for a stationary 
surface can therefore be written as 



It can be seen from equation (l8) that the average net heat-transfer 
coefficient for the oscillating surface would be larger than that for 
the stationary surface (that is, Rg > l) only when T W >T„, and it would 
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fee lower when T^< T„. For a given temperature difference, therefore, 
the ratio of the average net heat-transfer coefficients, as given by- 
equation (l8), can be altered by changing the amplitude of the 
oscillations . 


If the primary interest is not in the net heat transferred across 
the surface (as it would be, for example, for regenerators) but is in 
the thermal state of the fluid itself, the effect of the oscillations 
qan be studied by comparing the average static temperature defect de- 
fined as 




(19) 


for the oscillating surface to that for the stationary surface. The 
ratio of these defects is given by 





3 3" 

^osc U^Pr 3^%^ . 

Prit 

.(*0 * ? f - 

^"sta ® nc p(!«“®v) | 

2n 


For the limiting 


case, tg -*■ 0, equation (19) reduces to 



( 20 ) 

( 21 ) 


It can thus be seen that males s Fr>64/9 the oscillations will not be 
advantageous for lowering the fluid temperature for the case when 


T. 


'V 


A more meaningful parameter to examine would perhaps be the ratio 
of the average total -temperature defect with the oscillations to that 
for conduction to the stationary surface . The average total -temperature 
defect is defined as in equation (l9) except that T is replaced by the 
total temperature T^, which is given by 
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In this way the thermal condition of the fluid independent of its motion 
is indicated. This ratio is given by 


Ra * 1 - 


U 2 Pr 


4c p (T. - T w ) J 


r 

2 n 





Thus, the ratio *4 could be greater than unity dependent on the rela- 
tive orders of magnitude of the terms in the braces in equation ( 22 ) 
and in the limiting case as tQ -*■ 0 . 



U^Pr 

3 (1 - Pr) 


- 8 -V s J 


(23) 


Thus, in this ease if T„>T W the fluid would be at a lower relative 
temperature with the oscillations (R 40 > l) only if Pr>9. 

An analogous discussion, of course, follows for increasing the 
fluid temperature. 


COECLUBUSG REMARKS 

The results of this preliminary analysis of the temperature distri- 
bution over conducting oscillating surfaces have shown that the heat 
transfer associated with the oscillating surface can be significantly 
different from that for conduction between the same initial temperature 
difference and a stationary surface. In addition to the heat transfer, 
the thermal state of the fluid was studied by means of average static- 
or total-temperature defects and it was found that the surface oscil- 
lations can alter the thermal state of the fluid. In each case, the 
range of parametric values for which the surface oscillations are bene- 
ficial were indicated. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 19, 1954 
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APPENDIX - SYMBOLS 

The following symbols are used In this report: 
specific heat at constant pressure 


B 


10 


R, 


R 


30 


R 

T 

T, 


'40 


w 


T 

A w 

T 

•0 

ZSD 


2 


S S 


fj 


exp 


(-r 2 )dr 


thermal conductivity coefficient 
frequency of oscillations 
Prandtl number 

ratio of average total heat-transfer coefficients 

ratio of average total heat-transfer coefficients evaluated as 

to -° 

ratio of average net heat-transfer coefficients 

ratio of average static-temperature defects 

ratio of average static-temperature defects evaluated as 
t Q -0 

ratio of average total-temperature defects 

ratio of average total-temperature defects evaluated as 

t 0 -° 

static temperature 

equilibrium (or steady-state) static temperature 
total temperature 
surface temperature 

surface temperature for steady state 
static temperature in ambient fluid 
temperature defect 
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t time 

tQ time at the start of an arbitrary cycle 

U amp litude of oscillations 

u velocity parallel to surface 

y coordinate normal to surface 

a thermal, diffusivity 

1 ) von Mises coordinate normal to surface 

p absolute viscosity 

v kinematic viscosity 

p density 

Subscripts: 

c complementary solution 

osc oscillating surface 

p particular solution 

sta stationary surface 
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